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Abstract 

The enhancement of the transitions Vn (v e ) ~^ v e {ya\r)t v 2 ~ > v ei &u, ~ > v s , 
v e — ► u s , etc. of neutrinos passing through a medium of nonperiodic density 
distribution, consisting of i) two layers of different constant densities (e.g., Earth 
core and mantle), or ii) three layers of constant density with the first and the 
third layers having identical densities and widths which differ from those of the 
second layer (e.g., mantle-core-mantle in the Earth), is studied. We find that, 
in addition to the known local maxima associated with the MSW effect and 
the neutrino oscillation length resonance (NOLR), the relevant two- neutrino 
transition probabilities have new resonance-like absolute maxima. The latter 
correspond to a new effect of total neutrino conversion and are absolute maxima 
in any independent variable: the neutrino energy, the layers widths, etc. We 
find the complete set of such maxima. We prove that the NOLR and the new 
absolute maxima are caused by a maximal constructive interference between 
the amplitudes of the neutrino transitions in the different density layers. We 
show that the strong resonance-like enhancement of the transitions in the Earth 
of the Earth-core-crossing solar and atmospheric neutrinos is due to the new 
effect of total neutrino conversion. 



1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 
1784 Sofia, Bulgaria 



1 Introduction 



It was pointed out in 0] that the v 2 — > u e , — > u e (u e — > i^( r )) oscillations of solar 
and atmospheric neutrinos in the Earth, caused by neutrino mixing (with nonzero 
mass neutrinos) in vacuum Q, can be strongly amplified by a new type of resonance- 
like mechanism which differs from the MSW one and takes place when the neutrinos 
traverse the Earth core on the way to the detector (see also 0, [|, ffl). As numerical 
calculations have shown, at small mixing angles (sin 2 29 ^ 0.05), the maxima in the 
neutrino energy E of the corresponding transition probabilities, P(yi — > v e ) = Pe2 
and P(Vfi — > v e ) {P(ye u (j.m))i caused by the new enhancing mechanism, are 
absolute maxima and dominate in P e2 and P[y^ — > u e ) (P(v e u tiM))' the values 
of the probabilities at these maxima in the simplest case of two-neutrino mixing are 
considerably larger - by a factor of ~ (2.5 — 4.0) (~ (3.0 — 7.0)), than the values of 
P e 2 and P(v ll — > Ve) = P[y e u h(t)) a t the local maxima associated with the MSW 
effect taking place in the Earth core (mantle). The effect of the new enhancement is 
equally dramatic at large mixing angles. The enhancement is present and dominates 
also in the v 2 — > v e transitions in the case of v e — v s mixing and in the v e — *> v s and 
v n —* v s transitions at small mixing angles 0, |3|, ffl. Even at small mixing angles 
the enhancement is relatively wide |2|, in the Nadir angle 0, h, and in the neutrino 
energy and leads to important observable effects. It also exhibits rather strong energy 
dependence. 

The new resonance-like enhancement of the probability P e i has important impli- 
cations for the tests of the MSW v e — > v^r) and u e — > v s transition solutions of 
the solar neutrino problem Jl], ||], |^, ||, 0, §]. For values of Am 2 from the small 
mixing angle (SMA) MSW solution region and the geographical latitudes at which 
the Super-Kamiokande, SNO and ICARUS detectors are located, the enhancement 
takes place in the u e — > u^r T \ case for values of the 8 B neutrino energy lying in the 
interval ~ (5 — 12) MeV to which these detectors are sensitive. Accordingly, at small 
mixing angles the new resonance is predicted ]5|, ^ to produce a much bigger - by a 
factor of up to ~ 6 (~ 8), day-night, or D-N, asymmetry in the Super-Kamiokande 
(SNO) sample of solar neutrino (charged current) events, whose night fraction is due 
to the core-crossing neutrinos (Core D-N asymmetry), in comparison with the asym- 
metry determined by using the whole night event sample (Night asymmetry). As a 
consequence, it can be possible to test a substantial part of the MSW v e — > io r \ 
SMA solution region in the Am 2 — sin 2 Id plane (see, e.g., flQfl ) by performing core 



D-N asymmetry measurements [^, 0, [| . The Super-Kamiokande collaboration has 
already successfully applied this approach to the analysis of their solar neutrino data 
|I0] : the limit the collaboration has obtained on the D-N asymmetry utilizing only 
the measured event rate in the night N5 bin, approximately 80% of which is due to the 
Earth-core-crossing solar neutrinos, permitted to exclude a small part of the MSW 
SMA solution region. In contrast, the published Super-Kamiokande upper limit on 

2 The V2 — > v e transition probability accounts, as is well-known, for the Earth effect in the solar 
neutrino survival probability in the case of the MSW two-neutrino v e — + v ^(t) an d v e — > v s transition 
solutions of the solar neutrino problem, v s being a sterile neutrino. 

3 The Nadir angle determines uniquely the neutrino trajectory through the Earth. 
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the Night D-N asymmetry [IIJ does not allow to probe the SMA solution region: the 
predicted asymmetry is too small (see, e.g., ||). 

The same new enhancement mechanism can and should be operative also [l] in 
the —>■ p e [y e — > z^t( T )) small mixing angle transitions of the atmospheric neutrinos 
crossing the Earth core if the atmospheric and indeed take part in large mixing 
vacuum v^ip^) <-> u T (u T ), oscillations with Am 2 ~ (10 -3 — 8x 10 -3 ) eV 2 , as is strongly 
suggested by the Super-Kamiokande data [H|, and if all three flavour neutrinos are 
mixed in vacuum. The new resonance can take place practically for all neutrino 
trajectories through the core (e.g., , for the trajectories with h = (0° — 23°)). It can 
produce an excess of e-like events at —1 < cos9 z < —0.8, 9 Z being the Zenith angle, in 
the Super-Kamiokande multi-GeV atmospheric neutrino data and can be responsible 
for at least part of the strong Zenith angle dependence of the Super-Kamiokande 
multi-GeV and sub-GeV \i— like data ||, ||. 

The Earth enhancement of the two-neutrino transitions of interest of the solar and 
atmospheric neutrinos at relatively small mixing angles has been discussed rather 
extensively, see, e.g., refs. |7|, [T2], [TJ, [TJJ]. Some of the articles contain plots of 
the probabilities P e 2 and/or Piv^ — > u e ) or P[y e — > v^( T )) on which one can clearly 
recognize now the dominating maximum due to the new enhancement effect. However, 
this maximum was invariably interpreted to be due to the MSW effect in the Earth 
core (see, e.g., [12|, [13|]) before the appearance of Q. 



In view of the important role the new mechanism of enhancement of the neutrino 
transitions in the Earth can play in the interpretation of the results of the experiments 
with solar and atmospheric neutrinos and in obtaining information about possible 
small mixing angle [v^] and v e (z/ e ) oscillations, it would be desirable to have an 
unambiguous understanding of the underlying physics of the mechanism. In |T| it was 
interpreted as an effect of constructive interference between the various probability 
amplitudes (notably of the neutrino transitions in the Earth mantle and in the Earth 
core), entering into the sum representing the probability amplitude of the transition 
in the Earth. 

The exact conditions for the enhancement of the probabilities P e 2, P{v^ — > v e ) 
(P(v e ^(r))' e tc by the new mechanism can be studied analytically [p]] owing to 
the results of a detailed numerical analysis |]15] (see also which showed that for 
the calculation of the probabilities of interest the two-layer model of the Earth density 
distribution provides a very good (in many cases - excellent) approximation to the 
more complicated density distributions predicted by the existing models of the Earth 
|16| , pTT[1 . In the two-layer model, the neutrinos crossing the Earth mantle and the 
core traverse effectively two layers with constant but different densities, p man and p c , 
and chemical composition or electron fraction numbers, YJ nan and Y" e c . The neutrino 
transitions of interest [] result from two-neutrino oscillations taking place i) first in 
the mantle over a distance X' with a mixing angle 9' m and oscillation length L man , 
ii) then in the core over a distance X" with mixing angle 9!^ ^ 9' m and oscillation 
length L c ^ L man , and iii) again in the mantle over a distance X' with 9' m and 
Lman- Utilizing the above prescription, analytic expressions for the probabilities of 



4 The densities p m an,c should be considered as mean densities along the neutrino trajectories. 
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neutrino transitions in the Earth P e 2, Piy^ — *> u e ) = P(u e — > ^(r)) ; etc. were derived 
and were used to study their extrema Q. Assuming that the neutrino oscillation 
parameters Am 2 > and cos 29 > are fixed and treating the neutrino paths X' 
and X" as independent variables, it was found in [|lj that in addition to the maxima 
corresponding to the MSW effect in the Earth mantle or in the Earth core, there exists 
a new maximum in P e2 and Piy^ —>■ u e ). The latter takes place when the relative 
phases the states of the two energy-eigenstate neutrinos acquire after the neutrinos 
have crossed the mantle, AE'X' = 2irX' /L man , and the core, AE"X" = 2nX"/L c , 
obey the constraints, 

AE'X' = ix(2k + 1), AE"X" = n(2k' + 1), k, k' = 0, 1, 2, ... , (1) 

and if the inequality 0, |3| 

cos(2C-A9' m + e) <0, (2) 

is fulfilled. Condition (2) is valid for the probability P e 2- When equalities (1) hold, 
(2) ensures that P e 2 has a maximum. At the maximum P e 2 takes the form Q 

P2 ax = ^\2d" m -Ae' m + e). (3) 

The analogs of eqs. (1) - (3) for the probability P{v^ e ) — > ^(^r)) can be obtained 
H by formally setting 6 = while keeping ff m ^ and 0'^ ^ in (1) - (3) Q. The 
term "neutrino oscillation length resonance" (NOLR) was used in to denote the 
resonance-type enhancement of P e 2, Piy^ — > Ve) {P{ v e ^(r)))j etc associated with 
the conditions (1) - (2). 

The new mechanism of strong enhancement of the probabilities P e 2, Piy^ — ► v e ) 
(P{y e v h{t)))i e tc. for the Earth core crossing neutrinos, was identified in U 
with the neutrino oscillation length resonance. This interpretation was based on the 
observation that for the several fixed test values of sin 2 26 ^ 0.02 and Nadir angle 
h considered in Q, conditions (1) were approximately satisfied at the corresponding 
non-MSW absolute maxima of P e 2 and Piy^ — > u e ) in the variable E/Am 2 and the 
numerically calculated values of P e 2 and Piy^ — * v e ) at these maxima were reproduced 
by eq. (3) with a relatively good precision. Doubts about the correctness of such an 
interpretation remained since i) for the parameters of the Earth [|16j (core radius, 
mantle width, etc.), used in the calculations, the set of conditions (1) - (2), as was 
shown in cannot be exactly satisfied at small mixing angles, and ii) although the 
phase AE'X' was close to 7r at the relevant maxima of P e 2 and Piy^ —* u e ) in the 
test cases studied, the values of the phase AE"X" differed quite substantially from 
those required by eq. (1). Moreover, the same (or similar) enhancement mechanism 
was found to be operative in the — > v e transitions in the case of v e — v s mixing and 
in the v e v s transitions, at small mixing angles However, as was noticed in [|TJ, 
the conditions (1) - (2) are not even approximately fulfilled for the z/ 2 — > v e = v s —>■ v e 

5 The conditions for the corresponding maximum in P{v e —> v s ) and P(v^ — > v s ) coincide in form 
with those for P(^( e ) -> v e ( K T))- 
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and v e — > z/ s transitions at small mixing angles; similar conclusions were reached for 
the — > i> s transitions ||. 

The indicated results stimulated the systematic study of the extrema of the prob- 
abilities P e 2, P{ffi — > ^ e ), f (^e ~^ ^s), etc. of the neutrino transitions in the Earth 
and more generally, in a medium of non-periodic constant density layers, which is 
presented here. We consider transitions of neutrinos in a medium consisting of i) 
two layers of different constant densities, and ii) three layers of constant density with 
the first and the third layers having identical densities and widths which differ from 
those of the second layer. Neutrinos pass through such systems of layers on the way 
to the neutrino detectors, e.g., when they i) are produced in the central region of 
the Earth and traverse both the Earth core and mantle, ii) travel first in vacuum 
and then in the Earth mantle, iii) cross the Earth mantle, the core and the man- 
tle again. We derive and analyze the complete set of extrema of the two-neutrino 
transition probabilities in both cases of a medium with two and three constant den- 
sity layers, assuming that U Am 2 /E and sin 29 are fixed and treating the widths 
of the layers, X' and X", as independent variables. For both media considered we 
find that in addition to the local maxima associated with the MSW effect and the 
NOLR, there exist absolute maxima caused by a new effect of enhancement and cor- 
responding to a total neutrino conversion. We find the complete set of such absolute 
maxima. The latter are absolute maxima in any possible independent variables: the 
neutrino energy, the widths of one of the layers, etc. The conditions for existence 
of the new maxima are derived and it is shown that they are fulfilled, in particular, 
for the transitions of the Earth-core-crossing solar and atmospheric neutrinos. These 
conditions differ from the conditions of enhancement of the probability of transitions 
of neutrinos propagating in a medium with periodically varying density, discussed in 



T8| . We show that the strong resonance- like enhancement of these transitions is due 
to the effect of total neutrino conversion. At small mixing angles and in the case 
of the transitions vi —>■ v e and — * u e , for instance, the values of the parameters 
sin 2 29 and Am 2 /E at which the total neutrino conversion takes place for neutrinos 
traversing the Earth core along a given trajectory are rather close to the values of 
the parameters for which the NOLR giving P e2 = 1 or Piv^ — > i/ e ) — 1 occurs, and 
the two enhancement mechanisms practically coincide. For smaller sin 2 29 the NOLR 
reproduces approximately (in some cases - roughly) the enhancement: for fixed h and 
sin 2 29 = (10~ 3 — 10 -2 ), for instance, it reproduces the values of P{yi — > v e ) and 
P{ u fi ~ * v e) & t the absolute maxima in the Am 2 /E variable with an error of ~ (15 - 
60)% (see further). In the case of the v% — > u e = u s — > u e , u e — > v s and the — > v s 
transitions at small mixing angles, the total neutrino conversion occurs for values of 
the parameters for which the NOLR is not realized. In all transitions, however, only 
the total neutrino conversion mechanism is operative. This mechanism is responsible 
also for the strong enhancement of the above transitions at large mixing angles as 
well as of the z/ 2 — * v e transitions due to v e — v s mixing and of the u e — > v s and the 
~~ > v s transitions at small and large mixing angles. 

We demonstrate that the NOLR and the newly found resonance- like enhancement 
effect are caused by a maximal constructive interference between the amplitudes of 
the neutrino transitions in the different density layers. Thus, the local and abso- 



4 



lute maxima they produce in the neutrino transition probabilities are of interference 
nature. Therefore we confirm the physical interpretation of the strong resonance- 
like enhancement of the transitions in the Earth of the Earth-core-crossing solar and 
atmospheric neutrinos, given in |]. 

Most of the results of our analysis are illustrated by the physically realistic exam- 
ples of the transitions of neutrinos passing through the Earth. These examples and 
the results regarding the transitions of Earth-core-crossing solar and atmospheric neu- 
trinos are derived using the Stacey model from 1977 |16|] as a reference Earth model. 
The density distribution in the Stacey model is spherically symmetric and in addition 
to the two major density structures - the core and the mantle, there are a certain 
number of substructures (shells or layers). The core has a radius R c = 3485.7 km, 
the Earth mantle depth is approximately R ma n — 2885.3 km, and the Earth radius 
in the Stacey model is Rq = 6371 km. Therefore, when the Nadir angle h is less 
than 33.17°, neutrinos arriving at the detector pass through three layers: the mantle, 
the core and the mantle again. The mean values of the matter densities of the core 
and of the mantle read, respectively: p c = 11.5 g/cm 3 and p ma „ = 4.5 g/cm 3 . The 
density distribution in the 1977 Stacey model practically coincides with that in the 



more recent PREM model [|T7 1 . For Y e we have used the standard values | 16| , [17], [19 1 



see also M) Y™ an = 0.49 and K c = 0.467. 



A brief description of the results of the present study is given in ref. f2TJ . 



2 Preliminary Remarks 

We will consider the simple case of mixing of two weak-eigenstate neutrinos v a and vp, 
a 7^ (3 = e, /i, r, s, v s being a sterile neutrino, in vacuum. The mixing matrix relating 
the states of the neutrinos v a ^p and of the neutrinos having definite masses 
in vacuum is chosen in the form: 



v a \ _ f cos 9 sin 9 \ f v\ 
up J I — sin 8 cos 9 J \ u 2 



(4) 



where 9 is the vacuum mixing angle. In the case of relativistic and stable neutrinos 
z/i 5 2, the evolution equation for neutrino propagation in constant-density medium 
reads 

. d_ / v a \ _ AE ( -cos(2# m ) sin(2# m ) \(v 0l \_ hA (v a \ 
dt\vp ) 2 ^ sin(2# m ) cos(2£ m ) ) [ v p ) 1 } 

Here 9 m is the mixing angle in matter, 

cos(20 m ) = -i- ((Am 2 /2£) cos(2^) - V a p), (6) 

E being the neutrino energy, AE is the difference between the energies of the two 
energy-eigenstate neutrinos in matter, 



A£=^(co S (2 e) -^) 2 W(2e), (7) 
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Am 2 = m\ — m\, and V a p is the difference between the effective potentials of v a and 
vp in the medium. We shall always assume in what follows that 

cos(2#) > 0, Am 2 > 0. (8) 

In the case of neutrino propagation in an electrically neutral unpolarized cold medium, 
like the Earth, one has: 

V e , = V2G F N e , V es = V2G F (N e -^N n ), V, s = -V2G F N n /2, (9) 

where N e and N n are the electron and the neutron number densities in the medium, 
respectively. The antineutrino states v a and Vp satisfy the same equations (5) with 
V a/ 3 replaced by V & p = —V a p in eqs. (6) and (7). For the Earth N e = N n and in 
addition to V efl > we have V es > and > 0. If (8) holds, the neutrino mixing can 
be enhanced by the Earth matter only in the cases of (v e ) — > v e (v^) (y 2 — > u e ), 
v e - > v s and — > v s transitions, while for, e.g., cos(26 l ) > and Am 2 < this will 
be true for the P M (P e ) — > v e {v KT ) [p2 — > v e ), v e — > v s and — > v s transitions. Our 
general results will be formulated for the generic weak-eigenstate neutrino transitions, 
v a vp\ for neutrinos crossing the Earth we will be interested either in the former or 
in the latter set of transitions, depending on whether eq. (8) holds or Am 2 cos(2#) < 0. 
The case of (solar) z/ 2 — > v e transitions in a three-layer medium (the Earth) will be 
considered separately. 

The evolution of the neutrino system is given by the unitary matrix (see, e.g., 
0,0) 

U = exp(—ij\4t) = cos0 — i(trn) sin0, (10) 

where 

0=±AEt (11) 

and 

n= (sin(2fl m ),0,-cos(2# m )) (12) 

is a real unit vector. The probabilities of the u a — > vp [yp — > u a ) transition, P a p(p a ) = 
\A a p(p a ) | 2 , and of the u a (up) survival, P aa (pp) = \A aa (pp)\ 2 , are determined by the 
elements of the evolution matrix U, which coincide with the four different probability 
amplitudes A a p = Up a , etc.: 

P a fi = \Up a \ 2 = (nisin0) 2 + (n 2 sin0) 2 = Pp a (13) 

and 

P aa = \U aa \ 2 = cos 2 + (n 3 sin 0) 2 = 1 - P a p = Ppp. (14) 

We are interested in the extrema of the transition probability P a p when neutri- 
nos propagate through a medium with (nonperiodic sequence of) layers of different 
constant density and chemical composition. In this case the evolution matrix, as is 
well-known, represents a product of the evolution matrices for the different layers and 
can always be written in the same form as in eq. (|TU]) . It follows from ( ]XB| ) and ( |HD 
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that the conditions for an absolute maximum of the transition probability, P a p = 1, 
read 

max P a(i = 1 : < C ° S ^ ° (15) 
' I n 3 sin = 0, 

while the conditions for an absolute minimum of P a p have the form: 

minP a/3 = 0: \ niSm ^ ^ (16) 
M I n 2 sin = 0. 

When neutrinos propagate in a constant-density medium the transition probability 
P a( i = \{l - cos 20) sin 2 (2# m ) = (sin0sin(2fl m )) 2 (17) 

can be non-negligible even in the case of small vacuum mixing angle 9 if the neutrino 
mixing in matter is maximal, 9 m = it/ '4. This is the well-known MSW effect which 
takes place when the resonance condition 

cos(2fl m ) = (18) 

is fulfilled. In order to get total neutrino conversion, P a p = 1, the additional require- 
ment 

cos0 = O (19) 

has to be satisfied. For a fixed time of propagation t, or a distance X = t traveled 
by the neutrino, there always exists a solution of the system of equations (p~5|) The 
absolute minima P a p = correspond, as it follows from ([if]), to the curves 

sin0 = O. (20) 

In the next sections we discuss the case of neutrino propagation through a number 
of nonperiodic const ant- density layers. A new phenomenon of maximal constructive 
interference between the probability amplitudes of the transitions in the different 
layers takes place in this case, leading to a substantial enhancement of the transition 
probability. Even when the oscillation parameters have values very different from 



those required by the resonance conditions (JlSf) for the individual layers, the neutrino 



transition probability can reach its absolute maximum, P a p = 1, due to this effect. 



3 Medium with Two Constant - Density Layers 

Suppose neutrinos v a or vp have crossed a medium consisting of two layers having 
constant but different density and chemical composition. These could be the Earth 
core and mantle (for neutrinos born in the Earth central region), the vacuum and the 
Earth mantle, etc. Let us denote the effective potential differences and the lengths of 
the paths of the neutrinos (antineutrinos) in the two layers by V' a p (Vis), O^g) 
and X', X", respectively. We shall assume without loss of generality that 

o< K pm \ < K^\. (2i) 
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The transitions of v a and vp crossing the two layers will be determined by the two sets 
of two parameters - the mixing angle in matter (6) or the unit vector (12), and the 
phase difference (11) with t — X' or X", which characterize the transitions in each 
of the layers: 6' m or n' and <fr', and 6'^ or n" and 0". Since the evolution matrix of 
the system is given by U = U"U', U' and U" being the evolution matrices in the first 
and in the second layers, the parameters of U, <£> and n (eq. (|10D), can be expressed 
in terms of (ft', n' and 0", n" as follows (see, e.g., [^lj ): 

cos $ = cos 0' cos 0" — (n' • n") sin 0' sin 0" 

= cos 0' cos 0" - cos(2C - 29' m ) sin 0' sin 0", (22) 
n sin $ = n' sin 0' cos 0" + n" cos 0' sin 0" — [n' x n"] sin 0' sin 0". 

Using, e.g., eqs. fllB]) and ( |22] ) it is not difficult to derive the transition probability 

Pa/3 — P/3a in this case: 

P a/3 = - [1 - cos 20'] sin 2 26' m + - [1 - cos 20"] sin 2 2Q" m + - [1 - cos 20'] [1 - cos 20"] x 

sin 2 (2C-2^)-sin 2 2^-sin 2 2Cj + \^2(j)' sin 20" sin 26' m sin 2Q" m . (23) 

Due to T-invariance and unitarity, the probability (^) is symmetric with respect to 
the interchange of the parameters of the two layers, i.e., does not depend on the order 
in which the neutrinos traverse them. 

For fixed sin 2 (2#) and Am 2 /E and given density and chemical composition in the 
two layers, the unit vectors n' and n" are also fixed for each layer. The phases 0' and 
0", which depend on the widths of the layers X' and X", can be independent variables 
of the system if the two widths X' and X" can be varied independently. This is an 
interesting case with a clear physical meaning and we are going to investigate it in 
the present and the next Sections. In this way the NOLR resonance in the v-i — > v e , 
v ti (^e) ~" *■ v e {. u h;t) and v e — > v s transitions was discovered |l|. Let us note that in 
the case of the Earth, X' and X" cannot be treated as independent variables as long 
as the Earth radius and the Earth core radius are fixed. Neutrino transitions in the 
Earth will be considered at the end of this Section and in Section 6. 

Varying the phases 0' and 0" we can investigate the number and the structure of 
the extrema of the transition probability P a p. These are determined by the system 
of two equations 



dPa/3 



dPct/3 



sin 26' m F(26' m - 26^ 2B" m , 20', 20") = 0, (24) 
sin 2B" m F{2Q" m - 29' m) 26' m , 20", 20') = 0, (25) 



where 

F(Y, Z; ip, ip) = sin Y cos Z sin (p + sin Z [sin ip cos ip cos Y + cos ip sin ip] , (26) 

and by two known supplementary conditions on the values of the second derivatives 
of P a p at the points where eqs. (f2~4| ) and ( [25] ) hold. 
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We are interested first of all in the absolute maxima of the neutrino conversion, 

case A : max P a/ 3 = 1. (27) 

They are determined by the equations flT5|) and ( P2"|) 

f cos $ = cos 0' cos 0" - cos(2C - 26 J sin 0' sin 0" = 
a/3 \ n 3 sin $ = - [cos(2%) sin 0' cos 0" + cos(20£) cos 0' sin 0"] = 0. 

(28) 

It is not difficult to check that conditions (|24]) and (p5|) are satisfied if relations (f28| 
hold. The solutions of (28) can be readily found: 



solution A : < 



tan0' = ±^ 
tan0" = ±- 



; - cos(2Q 

cos(2^) 008(2^^2^1 ' (29) 

I -cos(2Q 1 ' 

' cos(2C) cos(2C - 26 J ' 



where the signs are correlated. Obviously, the solutions ( p9|) do not exist for any pairs 
of values of the parameters sin 2 (2#) and Am 2 /E. Under the assumptions (§), (pl|) 
and for V^g > (e.g., for the z/^ (z/ e ) — > v e (v^ T ), v e — > v s and — > v s transitions of 
the Earth-crossing neutrinos), they can take place only in the region limited by the 
three conditions (Fig. 1) 

f cos2^>0 

region A: I cos 26'^ < (30) 
{ 008(2% - 26'J > 0. 

The first boundary, cos 26' m = 0, corresponds to a total neutrino conversion in 
the first layer. From ( p9{) we get the correct result for the relevant conditions on the 
phases (see (|19|) and (pOl)), 

, ,. R / cos 0' = O, or 20' = 7r(2fc' + l), fc' = 0,1,..., , . 

solution B : ( . v , / (31) 

I sm0 = 0, or 20 = 27TK , re" = 0, 1, 

guarantying that i) neutrino conversion takes place in the first layer, and ii) no neu- 
trino conversion occurs in the second layer, so that when cos 2d' m = we have maximal 
transition probability, P a p = 1. It is easy to show by direct calculations using eqs. 
that the phase conditions ( Plf ) correspond to the local maxima 



case B : max P a p = sin 2 26 J (32) 
provided the oscillation parameters sin 2 (2#) and Am 2 /E belong to the region B, 

region B : cos2%<0. (33) 

6 These conditions can also be derived from the two supplementary inequalities ensuring that the 
solutions (29) of (EJ) and rfgq) indeed correspond to maxima of P a p- 
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It proves useful to interpret this result in terms of probability amplitudes. If 
neutrinos v a have crossed the two layers, the amplitude of the v a — » v$ transition 
represents a sum of two terms: 

A aP = U Pa = U'> a U' aa + V%U'^ (34) 

where 



U' aa = cos <P' + i cos(20 sin 0', = -i mn(20£) sin 0", 

= -i sin(2^) sin 0', Ufa = cos 0" - i cos(20^) sin , 



(35) 



In the general case the transition probability 

Pop = W\U' aa \ 2 + WW + 2Re ((u^U' aa )* U^U' pa ) (36) 

is a sum of two products of the probabilities of neutrino oscillations in the different 
layers and of the interference term. The latter can play a crucial role in the resonance- 
like enhancement of the neutrino transitions. 

The phase requirements ( |3"1~D reduce the expression in eq. fl3"6| ) to only one term, 



P a /3 = I^J 2 ) with no contribution from the interference term because Up a = 0. Thus, 
the maxima in region B (|33|) can be ascribed to neutrino transitions taking place in 
the first layer only. The absolute maxima, max P a p = 1, correspond for sin 2 29 < 1 
to the MSW effect in the first layer. 

The second boundary cos(26^J = in ( p0|) corresponds to a total neutrino con- 
version in the second layer, 

, .. n \ sin 0' = O, or 20' = 2ick', k' = 0,1,..., , Q7 s 
solution G : j = ^ ^ ^ = ^ + ^ = ^ ^ (37) 

This case is completely analogous to the previously considered case B, with the two 



layers interchanged. Using eqs. (|23| ) - ( p5|) we get solution C which realizes the local 
maxima 

case C : max P a p = sin 2 29'^, (38) 

in the region 

region C: cos20^>O. (39) 

The case of equality in eq. (|39|) can be associated with the MSW effect taking place 
in the second layer. 

We get a very different mechanism of enhancement of the neutrino transition 
probability P a p in the intermediate region 



cos 26' m > 0, 
cos 29" < 0, 



(40) 



located between the regions B ( p3|) and C fl39|). The maxima of P a p in this region 
are caused by a maximal contribution of the interference term in eq. (^) for P a p. 
Indeed, using eq. fl3~4] ) one can write the probability P a p in the form: 



a/3 



zi + z 2 | 2 , (41) 
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where Zl = (Re(U'^ a ), lm(U$ a U^)) and z 2 = (Re(U^U' 0a ), Im(^^J) are two 
vectors in the complex plane. In the case of the probability P a p the interference is 
maximal and constructive, as it follows from (41), when zi and z 2 are collinear and 
point in the same direction, i.e., when 



lm(U% a U' aa ) _ Re(U'^U' aa ) 



> 0. (42) 



The second equation in (|T5|) (or (^)) ensures that the two vectors are collinear, while 
the constraints ([5]) guarantee that they point in the same direction. Actually, the 
equation in ( f42"D coincides with the second equation in ([15]) for the resonance condition 

cos(2# m ) = -— [cos(20') sin <p' cos 0" + cos(2fl") cos <p' sin 0"] = 0, (43) 



with the additional constraints (f40|) , corresponding to the intermediate region where 
the new type of enhancement takes place. The transition probability of interest ( ^3|) 
can obviously be represented in the two layer case under discussion also in the form 

(0), 

P^ = sin 2 (2# m )sin 2 0, (44) 

where the first and the second multipliers are defined by eq. ([13|) and the first equation 
in (|^). To get a total neutrino conversion, P a/ 3 = 1, not only eq. (|43|) has to hold, 
but also condition ([19]) must be satisfied. This is possible in region A (|30|), but not 
in the whole intermediate region ([40|). 

Indeed, the boundary cos(26^ — 29' m ) = in (^) corresponds to a maximum of 
P a/ 3 associated with the phase constraints 

... J cos0' = 0, or 2(f)' = n(2k' + 1), k' = 0, 1, . . 

solution D : < „ „ \ „ .„ (45) 

I cos0 = 0, or 2<p = 7T{2k + 1), k = 0, 1, .... 

The above conditions are necessary for a maximal neutrino conversion in each layer 
(see (|i~9|) and (|3~5])). As can be shown exploiting eqs. (|23|) - (|^), they lead to local 
maxima 

case D : max P aP = sin 2 (2^ - 26' m ), (46) 
if the oscillation parameters belong to the finite region D, 

region D : cos(20^ - 26'J < 0. (47) 

Our results show a very interesting feature of the neutrino transitions in a medium 
consisting of two constant-density layers. One can have total neutrino conversion, 
P aj 3 = 1, even if the MSW resonance does not take place in any of the two layers. As 
we have seen, this is possible in region A, excluding the two boundaries cos(26^J = 
and cos(2#'^) = 0, corresponding to the MSW effect. Therefore, in order to have a 
large transition probability, P a p = 1, a periodic density profile is not required even 
when the MSW resonance is not realized in any of the two layers. These conclusions 
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remain valid, as we shall see, for transitions of neutrinos traversing three layers of con- 
stant density, e.g., for the transitions of the solar and atmospheric neutrinos crossing 
the Earth core on the way to the detectors. It should be emphasized that solution A, 
eq. (|29|) , providing the absolute maxima of P a p, can be and was derived directly from 
eqs. (15) and (22) without utilizing the standard extrema equations (24) - (25), i.e., 
its derivation does not depend on the chosen set of variables of P a/ 3, which are varied 
to obtain the extrema conditions. Obviously, it can also be obtained for a given set 
of variables from eq. (23) and the corresponding extrema conditions of the type of 
(24) - (25). " 

The absolute minima of P a p in the case under discussion are located on the inter- 
sections of the curves 

f sin0' = 0, or 20' = 2nk', k' = 0, 1, 
mm f a p - U : j ^„ = ^ or ^ = fc „ = ^ ^ (4Sj 

Our results apply to the case when the first layer is vacuum V = (Fig. 2), 
because there is no principal difference between the neutrino oscillations in constant- 
density medium and in vacuum. If we assume the nonzero density layer to be the 
Earth mantle, maximal neutrino conversion (solution A) can take place provided 
neutrinos travel a distance of ~ fewxlO 3 km in vacuum. 

The above results are valid also, e.g., for neutrinos born in the central region 
of the Earth. Such neutrinos cross two layers - the Earth core and mantle, on the 
way to the Earth surface. If the dimensions of the region of neutrino production 
are small compared to the Earth core radius, the neutrino path lengths X' and X" 
are fixed. The extrema condition with respect to the variable E is given in ref. 
It differs from the the extrema conditions (24) and (|25|). Accordingly, the 



solutions of the system (^J) - fl25|) do not correspond, in general, to extrema in the 
variable E. More specifically, solutions B, C and D, corresponding to local maxima 
of P Qj 3, max P a /3 < 1, in the variables and 0", do not provide local maxima 
in the variable E. Only solution A provides the absolute maxima of the neutrino 
transition probabilities of interest, P a p = 1, in any variable. Solution A is realized for 
v n { u e) — ► v e (^Vr); v e ~~ ► v s and — > v s transitions of neutrinos born in the central 
region of the Earth, as Table 1 demonstrates. 



4 Special Case of a Three - Layer Medium 

Similar analysis can be performed for the transitions of neutrinos traversing three 
layers of constant density and chemical composition when the first and the third 
layers have the same density, chemical composition and width which differ, however, 
from those of the second layer. This corresponds to the physically important case 
of solar and atmospheric neutrinos crossing the Earth core on the way to detectors 
located near or on the Earth surface, the first and third layers being the Earth mantle, 
and the Earth core playing the role of the second layer. 

We will use the same notations for the parameters of the first (third) and the 
second layers as in Section 2, and will assume that eq. ( pT|) holds. The evolution 
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matrix in the case of interest is given by: U = U'U"U' . Accordingly, the parameters 
of U, eq. (fL0|), can be expressed in terms of the parameters of the first (third) and 
the second layers as follows (see, e.g., f2lfl): 



f cos = cos(20') cos 0" - cos(20£ - 29' m ) sin(20') sin 0", , . 

\ nsin0 = n' [sin(20') cos0" - (n' • n")(l - cos(20')) sin0"] + n"sin0". ' ' 

An expression for the probability P a p is given in |jj (see also @). The necessary 
conditions for a maximum of the probability P a p, obtained by varying the phases 0' 
and 0", read 

f sin(2^) F(2C - 2^, vr/2; 0", 20' + tt/2) = 0, 
^ ' 1 F(JC - 29' m , 26> m , 0" + vr/2, 20') = 0, [M) 

while the requirements for an absolute maximum have the form 

mnrP f cos = F(2C-2C tt/2; 0", 20' + tt/2) = 0, , . 

Q/3 1 n 3 sin = F(2^-2C 2C + tt/2; 0", 20') = O. ^ 



where the function F(Y, Z;(p,ip) is defined by eq. (p6|). Using the expression for 
F(Y, Z; (p,ip) we get the conditions (|5T| ) in explicit form: 

p _ 1 J 2 cos(0') COS$ — COS0" = 0, s 
max ^ - : < cog ( 2 0//j + 2 sin ^/) cos ( 2 ^) cos $ = 0, 1 j 



where cos $ is given by eq. 

The absolute maxima of P a/ 3, 

case A : max P a p = 1, (53) 

are provided by the solutions of eq. ([y]) (or (|52|)), which can be found explicitly: 



solution A : < 



tan0' = ±4 



cos 20" 



v oo8(2^ - 4&J ' 

cos(2C) cos(2C - 40 ' 



where the signs are correlated. 

The probability P a p (Pan) exhibits a system of of maxima which is similar to that 
in the two layer case. Under the conditions (8), (21) and if V a a > (i.e., for the 

(z/ e ) — > v e (v KT ), v e — > v s and — > z/ s transitions in the Earth), solutions (|5"3] ) are 
realized in the region A (Figs. 3-5), 

reaion A ■ J cos ( 2 ^) ^ °' f55l 
region A . j _ ^ > Q (55) 

On the line belonging to region A, 

region B : cos 29' m = 0, (56) 
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we have 

case B : max P a p = sin 2 20' m = 1, (57) 

provided 

... / cos20' = O, or 20' =f(2fc'+l), A;' = 0,1,..., 

solution B : \ ,, „ 2 \„ „ (58) 

I sin =0, or 2<p — 2irk , k = 0,1, ... 

Besides these absolute maxima, there exist two regions, 

region C : cos(20£) > 0, (59) 



and 

with maxima 
and 



region D : cos(20£ - 40^) < 0, (60) 



case C : max P Qj g = sin 29'^, (61) 



case D : max P afi = sin 2 (20^ - 40^), (62) 
which correspond to the solutions 

, . . n J sin 0' = 0, or 20' = 2ttA;', k' = 0, 1, 
solution C: < „ „ , , „ s , „ (63) 

I cos0 = 0, or 20" = 7r(2/c" + 1), fc =0,1,..., 

and 

/ cos 0' = 0, or 20' = 71(2^' + 1), k' = 0, 1, 
somtlon ^ ' | COS0" = 0, or 20" = 7i{2k" + 1), jfc" = 0, 1, 1 > 

respectively. In contrast to the two-layer case, the region B is just a line, belonging 
actually to the region A, and on it only absolute maxima can be realized due to 
the MSW effect in the first (third) layer (the mantle in the case of the Earth). The 
case C corresponds to the MSW effect in the second layer, while the cases A and 
D correspond to the new resonance-like effect of constructive interference between 
transition amplitudes in the first (and third) and the second layers. Due to this effect 
the transition probability can reach its maximal value, P a p = 1, if the oscillation 
parameters sin 2 (20) and Am 2 /E belong to the region A, eq. (|55|). 

The absolute minima of the probability P a p are determined by the equation 

mm P a/3 = 0: F(26^ - 26' rn , 26' m , 0", 20') = 0, (65) 

while for F(29'^ n — 26' m , 26' m ; 0", 20') ^ 0, the necessary conditions for the minima are 
given by eq. (^D|) . 

We get the same system of maxima also in the v-i — > v e transition probability, 
P{i>2 —> v e ) = P e 2, V2 being the heavier of the two mass eigenstate neutrinos in vac- 
uum, which can be used to account for the the Earth matter effects in the transitions 
of solar neutrinos traversing the Earth. The probability P e 2 is given in the case of 
v e — io T ) mixing by the £7 e2 element, P e2 = \U e2 \ 2 , of the evolution matrix 
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U° = ( V" ) = ( V" V* ) ( C ° S ^ ^ ) ^ UO(9), (66) 

V U^i Up J \ U^ e U m J \ - sine cos6» J 

where U = U'U"U' . If v e — v s mixing takes place, the index /i appearing in the above 
equation has to be replaced by s. The matrix U & can be written in the general form 

my- 

U & = U'U'V'O = cos 0® - i(crn©) sin 0® . (67) 

It is not difficult to express the parameters of the evolution matrix (0), 0® and n®, 
in terms of the parameters 0', n', 0", n" and 9. The transition probability P e 2 can 
be written in the form 

P e2 = [F(29l - 29' m , 29> m - 9; 0", 20')] 2 + sin 2 9 \F{29" m - 29' m , vr/2; 0", 20' + vr/2)] 2 . 

(68) 

Using eqs. (|5"8"D and ([26]) we get the necessary conditions for the extrema of P e 2, 



max r aP , mm r aP . j F ^ _ ^ ^ _ Q _ ^ + ^ ^ = ^ 

while from eqs. ( |15j ) and (|16D , utilizing eqs. ([49"D and (|67D , one obtains the require- 
ments for absolute extrema of the probability P e2 : 

mnT P _ i . / cos <^ = ^( 2 C - 2 C vr/2; 0", 20' + vr/2) = 

maX ^ 2 " • \ nf sin 0° = F(2^ - 29> m , 29' m - 9 + vr/2; 0", 20') = 0, [ ( "> 

/ nf sin0® = F(2^ - 2^, tt/2; 0", 20' + tt/2) = 
mm 62 " 1 n§> sin 0® = F(2^ - 29' m , 29> m - B; 0", 20') =0. 



Conditions ([70]), for instance, have the following explicit form: 

J 2cos(0') COS$ - COS0" = 0, , - 

max ^2 - : j cos ( 2 ^ - 0) + 2 sin(0') cos(2^ - 9) cos $ = 0. (7 } 

These all are the same type of equations as (0) and we can easily find their solutions. 
The absolute maxima of P e 2, 

case A e : max P e 2 = 1, (73) 

are reached for 



solution A 



tan0' = ±* 



/ -cos(2<%-0) 
cos(2C-4^ + 6 

tan0" = + 



(74) 



cos(2C - 0) cos(2C - A9' m + 0) 
where the signs are correlated, in the region A & (Fig. 4), 

reazon A & ■ I C ° s(2 ^ " ^ " ° (75) 
region A . < _ + ^ > Q _ (75) 



15 



There are the line 

region B & : cos(2^ - 9) = 0, (76) 
belonging to A, and two bordering regions, 

region C & : cos(20£ - 0) > 0, (77) 

and 

re^on D° : cos(20^ - 40^ + 0) < 0, (78) 



where the solutions B, eq. (p8|) ; C, eq. (|63|) ; and .D, eq. ([64]), are realized. These 
solutions correspond to the following maxima of P e2 : 

case B & : max P e2 = sm 2 (29' m - 9) = 1, (79) 

case C® : max P e2 = sin 2 (20^ - 0), (80) 

and 

case D & : max P e2 = sin 2 (20^ - A9' m + 0). (81) 

Solutions D & and L> (eqs. ((0), (§2J), (0), QTH|) , (0)) correspond to the NOLR 
discussed [] in (l[|. 

For the absolute minima of P e2 we get the following solutions: 



tan (f)' = ±. 
tan 6" = ± 



sin (29 



sin(2C - 49' m + 9) 
sin(2C - 



(82) 



min P e2 = : < 

tan</>" = ±- 

^/sin^C - 0) sin(2C - < + 0) ' 

where again the signs are correlated. 

5 Transitions of Neutrinos Traversing the Earth 
Core 

In the case of transitions in the Earth of (solar and atmospheric) neutrinos which 
pass through the Earth mantle, the core and the mantle again, the two lengths X' and 
X" are not independent due to the spherical symmetry of the Earth, if both the Earth 
radius and the Earth core radius are fixed. We can choose as independent variables, 
for example, cos a, h being the Nadir angle, and y = Am 2 /E. The corresponding 
necessary conditions for the maxima of the probability P a p read: 

~ 0: IFF AE R®cosh -r-r-2(t) = 0, (83) 



dcos/i ' X" \ w d0" r d(20') 

dPap = . , y-2V^cos(29) dF X^_ y - 2V^cos(2fl) dF 

dy ' AAE' d(20') 2 AAE" d0" 



7 Solutions B and C (eqs. ( p8| ) and (|63|)) were considered briefly in pi as well. 
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V 

a 

AE 



+ sin(2#) sin 0' [cos(46' m - 29'^) sin 0' sin 0" - cos(2^) cos 0' cos . 



V'U sin(20) . 

smi 



XE , rl 2 ■ - [cos« - 2C) sin 2 0' + cos(2C) cos 2 0'j = 0, (84) 

where F = F{29'^ - 29' m , 29' m ; <f>", 20') is given by eq. ©, dF/d0" = F(2^ t - 
2C,2C;0" + vr/2,20') and dF/d(20') = sin2^(2C - 2^, vr/2; 0", 20' + tt/2). 
It is clear, that, in general, the maxima of P a p in the variables 0' and 0" do not 
correspond to maxima in the variables cos/i and Am 2 / E (compare eqs. fl50|) with 
eqs. (83) - (84)). It is not difficult to check, however, that solutions A, eq. ([54]) , 
for the absolute maxima, corresponding to a total neutrino conversion, P a p = 1, 
are solutions of the system of equations (83) - (84) as well. They give the absolute 
maxima of the neutrino transition probability P a p in any variable. Equation ( [54] ) 
gives also the complete set of such solutions. Note that solution A coincides on the 
curves cos(2^) = 0, cos(2^) = and cos(2^ - W m ) = with B, C and D, 
respectively. However, solutions C and D, eqs. (^) and fl64|), for the local maxima 
with max (Pap) < 1 no longer correspond to extrema in the variables cos h and 
Am 2 /E. New solutions for the local maxima of P a p, which do not coincide with 
the solutions associated with phases equal to multiples of n, are possible. The same 
conclusions are valid for the solutions A Q , C & and D & , eqs. ([74]), (|63| ) and ( |64] ) 
(giving (|30"D , (|8l]) ), for the absolute and local maxima of the probability P e 2- Our 
numerical studies confirm these conclusions. 

The solutions A, eq. (|54"D (A & , eq. (0)), providing a total neutrino conversion, 
Pa/3 — 1 {Pe2 — 1), exist for the Earth density profile and for all neutrino transitions 
of interest, v 2 — * z/ e , z/ M — > u e , v e — > i^( T ), ~ ^ and z/ M — > u s , and at small, 
intermediate and large mixing angles. They are realized for very different sets of 
values of the phases 20' and 20", which, are not multiples of tt. In Tables 2 - 4 we 
give a rather complete list of these solutions for the transitions u 2 — > v e , — > v e , 
v e v pl{t) 1 v e ~^ v s and — > v s of the Earth-core-crossing neutrinos, while Tables 5 - 
6 illustrate the behavior of the probabilities P^ and P e2 in the region of these solutions 
at small sin 2 29 (see further). The maximal neutrino conversion solutions A & and A, 
eqs. ([74]) and ([54]), are responsible for the strong resonance-like enhancement of the 
v 2 —* v e , i-V ~~ ► v ei v e u nM: v e v si etc. transitions in the Earth of the Earth- 
core crossing solar and atmospheric neutrinos, discussed in JT|, |2], |3|, f|]. In [JXJ] this 
enhancement was interpreted to be due to the neutrino oscillation length resonance 
- solutions D, eq. (|64]), for the neutrinos passing through the Earth core. At small 
mixing angles the values of the parameters at which the maximal neutrino conversion 
takes place for the v 2 — > u e — — > v e and v e — > transitions are rather close to the 
values of the parameters for which the NOLR giving P a p = 1 occurs (Tables 2 and 
4), while for the v 2 — > v e transitions caused by v e — v s mixing and the v e — > v s and 
Vfj, — > v s transitions these values are very different ]l|, [| (Tables 3-4). In both cases 
of transitions, however, only the maximal neutrino conversion mechanism is operative 
for the Earth-core-crossing neutrinos. 

Let us consider the case of v 2 — > v e , ^ —* v e [y e — > v^) transitions at sin 2 29 ^ 0.10 
in somewhat greater detail. As it follows from Table 2, for, e.g., h = 13° we have 
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P e2 = 1 at sin 2 26 = 0.044 and Am 2 /E = 7.1 x 10 -7 eV 2 /MeV, which corresponds to 
20' = 0.957T and 20" = 0.977T. At the indicated point one finds also that sin 2 (26^ — 
A6' m + 6) = 1, i.e., the NOLR solution D represents a very good approximation to 
solution A. Further, for any fixed sin 2 26 < 0.044, P e 2 has an absolute maximum in the 
variable Am 2 /E, which satisfies max P e 2 < 1 (see Fig. 7 and, e.g., Fig. 1 in |1|, Fig. 
4b in ||). When sin 2 26 increases from, say 10 -3 to 0.044, these absolute maxima form 
a continuous curve - a "ridge", with respect to the variables sin 2 26 and Am 2 /E, which 
passes through the total neutrino conversion point max P e2 = 1. It was suggested 
in I] on the basis of several test cases studied that the values of P e i forming the 
"ridge" are given with a relatively good precision, by the NOLR, i.e., by eq. (3) (or 
(jSl|)). These values have been calculated numerically and compared in Table 6 with 
the one predicted by eq. (3). We see from Table 6 that for sin 2 26 = 9.0 x 10~ 3 , eq. 
(3) reproduces the value of max P e i of interest with an error of 32%. For h = 23° 
and sin 2 26 = 10~ 2 the same error is 18%. The error increases with the decreasing of 
sin 2 26 and at sin 2 26 = 5.0 x 10~ 3 , for h = 13°; 23° reaches 43%; 23%. Results of a 
similar analysis for the probability P eil = P^ e are collected in Table 5. As it follows 
from Table 5, for h — 0°; 23° the NOLR, eq. (|64"D, describes the values of max P^ e 
on the corresponding "ridge" (see Fig. 6 and Figs. 5a - 5c in |Q) with an error which 
is 24%; 14% at sin 2 26 = 0.01 and increases to 42%; 19% at sin 2 26 = A x 10~ 3 . More 
generally, Tables 5 and 6 show that for fixed sin 2 26 = (10~ 3 — 10~ 2 ) and fixed h, 
the NOLR describes the absolute maxima of P e 2 and P^ e in the variable Am 2 /E, 
max P e 2(fie) < 1, with an error which varies between 14% and 60%. 

In some cases the total neutrino conversion solutions A and A° of interest occur 
for values of the parameters which are close to those for the solutions B, eq. (|58|), or 
C, eq. fl5"3"| ) (Tables 2-4). In all such cases the corresponding absolute maxima of 
the neutrino transition probabilities lie in the regions of the solutions A and A°, eqs. 
(H)and©. 

As Tables 2-4 and Figs. 6-9 indicate, for the Earth-core-crossing neutrinos, 
the new enhancement mechanism produces at sin 2 26 ^ 0.10 one relatively broad 
(in sin 2 26, Am 2 /E and h) resonance-like interference peak of total neutrino con- 
version in each of the probabilities P e 2, Piy^ — > v e )i P{v e ~~ > v s ) and Piv^ — > v s ). 
The total conversion peaks of the different probabilities are located in the inter- 
val 0.03 ^ sin 2 26 ^ 0.10. They all take place for values of the resonance density 
jyres _ j\ m 2 cos 26/ (2Ey2Gp), which differ from the number densities in the core 
and in the mantle V^p/ '(v2Gf) and V^ /3 /(\ / 2Gf), and lie between the latter two [0: 
V^p < \/2GFN res < V^p. For all the neutrino transitions considered, the correspond- 
ing interference peak is sufficiently wide in all variables, which makes the transitions 
observable in the region of the enhancement 0: the relative width of the peak in the 
neutrino energy is AE/E max = (0.3 — 0.5) and practically does not vary with sin 2 26 
(see also 0, Q), while in sin 2 26 it is ~ (2.0 — 2.3); the absolute width in h for the 
peaks located at h 25° is ~ (25 — 30)°. The points on the "ridges" leading to these 
peaks at, say, sin 2 26 ^ 0.03 represent for fixed h and sin 2 26 the dominating maxima 
in the variable Am 2 /E of the probabilities of interest, with max P a ^ < 1 (see Figs. 

8 Let us note that the peak is not symmetric in any of the variables. 
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6 - 9 and , e.g., Figs. 1 - 2 in and Figs. 4b, 5a - 5c, 10a - 10b, 15a - 15c in §). 

As Tables 2-4 and Figs. 6-9 show, the absolute maxima of total neutrino 
conversion, solution A, are present in the probabilities P e 2, P[y^ — > v e ), P{v e —> v s ) 
and Pipfj, — > v s ), at large values of sin 2 29 = (0.50 — 1.0) as well. They are also 
present in certain transitions at sin 2 29 = (0.15 — 0.50). Although these resonance- 
like interference maxima are, in general, narrower than at small mixing angles (see 
Figs. 6-9), they are sufficiently wide and can have observable effects both at sin 2 29 = 
(0.50 — 1.0) and sin 2 29 = (0.15 — 0.50). For certain specific trajectories the values 
of the phases 20' and 20" at the points of total conversion are close to odd multiples 
of 7r and the NOLR solution D reproduces approximately the results corresponding 
to the exact solution A. This takes place i) for P Me at h — 23° and sin 2 29 = 0.848 
(Table 2), ii) for P e2 (u e - z/ M mixing) at h = 23° and sin 2 29 = 0.941 (Table 4), 
iii) for P es at h = 0° and sin 2 29 = 0.999 and at h = 30° and sin 2 29 = 0.985 
(Table 3), iv) for P(z/ M -> v a ) at h = 30° and sin 2 29 = 0.981 (Table 3). This 
approximate equality of the two solutions is unstable with respect to the change of 
the parameters. Inspecting the Tables 2-4 one finds that for Nadir angles different 
from those indicated in i) - iv), the total neutrino conversion (solution A) takes place 
for values of the parameters at which the NOLR does not provide at all, or provides 
a poor description of the enhancement of interest. Note also that at sin 2 29 ~ 1 one 
can have sin 2 (26^ — 46^J = 1 at the points of total neutrino conversion, solution A, 
even for 20' and 20" which differ quite substantially from being odd multiples of 7r 
(as are the cases, e.g., of max P Me = 1 at h — 0° and sin 2 29 = 0.994, as well as of 
max P es = 1 at sin 2 29 = 0.931 and of max P^ = 1 at sin 2 29 = 0.992, for h = 23°). 
In these cases the analytic expression for P a/3 does not coincide with sin 2 (20^ —A9' m ). 

5.1 Transitions of Atmospheric Neutrinos 

— * v ei v e~^ u fj.(r)- These can be, e.g., sub-dominant transitions of the atmospheric 
Ufj,, driven by the values of Am 2 suggested by the Super-Kamiokande atmospheric 
neutrino data [ pT|j , 

Am 2 = (10~ 3 - 8 x 10" 3 ) eV 2 , (85) 

and by a relatively small mixing |TJ (see also 0, |4J). Such transitions should exist 
if three-flavour-neutrino (or four-neutrino) mixing takes place in vacuum, which is a 
very natural possibility in view of the present experimental evidences for oscillations of 
the flavour neutrinos. For the Earth-center-crossing neutrinos there are two solutions 
of the type A, providing a total neutrino conversion, P^e = Pep = 1, at small mixing 
angles (Table 2): at sin 2 29 = 0.034; 0.15 and Am 2 /P = 7.2; 4.8 x 10" 7 eV 2 /MeV. 
The first is reproduced approximately by the NOLR solution D while the second is a 
new type A solution. The above implies that for Am 2 = 10~ 3 eV 2 , the total neutrino 
conversion occurs at E = 1.4; 2.1 GeV, while if Am 2 = 5 x 10~ 3 eV 2 , it takes place at 
E = 7.0; 10.5 GeV. Thus, if the value of Am 2 lies in the region (85), the new effect of 
total neutrino conversion occurs for values of the energy E of the atmospheric v e and 
U/j, which contribute either to the sub-GeV or to the multi-GeV samples of e— like and 
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fi— like events in the Super-Kamiokande experiment. The implications are analogous 
to the one discussed in [p], |2|, U in connection with the NOLR interpretation of the 
enhancement of P^ e {P e p)- The new effect can produce an excess of e-like events in 
the region —1 < cos# z ^ — 0.8, 9 Z being the Zenith angle, in the multi-GeV (or a 
smaller one - in the sub-GeV) sample of atmospheric neutrino events, and should be 
responsible for at least part of the strong Zenith angle dependence, exhibited by the 
/i— like multi-GeV (sub-GeV) Super-Kamiokande data. 

The total neutrino conversion can take place at several values of Am 2 /E at large 
mixing angles (Table 2), sin 2 29 ^ 0.8. This can have implications, in particular, 
for the interpretation of the Super-Kamiokande data on the sub-GeV e-like events 
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— » v B . A total neutrino conversion due to the solution A, eq. (0), takes 
place at Am 2 /E < 10~ 6 eV 2 /MeV both at small and large mixing angles. For the 
Earth-center-crossing neutrinos the absolute maxima, Pp S = 1, are realized at (Table 
3) sin 2 26 = 0.07; 0.59; 0.78; 0.93; 0.999 for Am 2 /E = 4.1; 2.8; 9.7; 6.7; 3.4 x 
10~ 7 eV 2 /MeV. For the statistically preferred value of Am 2 = 4x 10" 3 eV 2 [|TTJ, 



this corresponds to E — 9.8; 14.3; 4.12; 6.0; 11.8 GeV, which is in the range of the 
multi-GeV fi— like events, studied by the Super-Kamiokande experiment. The total 
neutrino conversion maxima are present at all h ^ 30°, but at values of sin 2 29 and 
Am 2 /E, which vary with h. This variation (and the corresponding variation of 2(f)' 
and 20") is quite significant at large sin 2 29. 

We would like to add the following remark in connection with the enhancement of 
the transitions discussed. In the article by Q.Y. Liu and A. Yu. Smirnov, Nucl. Phys. 
B524 (1998) 505, it was noticed that in the case of muon (anti-) neutrinos crossing 
the Earth along the specific trajectory characterized by a Nadir angle h = 28.4°, 
and for sin 2 29 2* 1 and Am 2 / E ^ (1 - 2) x 10~ 4 eV 2 /GeV, the ^ -> v s transition 
probability is enhanced f\. The authors interpreted the enhancement as being due 
to the conditions 20' = tt and 20" = ir, which they claimed to be approximately 
satisfied and to produce dominating local maxima in the relevant transition probability 
P{vn v s ). Actually, for the values of the parameters of the examples chosen by 
Q.Y. Liu and A. Yu. Smirnov to illustrate this conclusion (Fig. 2 in their article) one 
has 20' = (0.6 — 0.9)7r and 20" = (1.2 — 1.5)7r. The indicated enhancement is due to 
the existence of a nearby total neutrino conversion point which for h = 28.4° is located 
at sin 2 29 S 0.94 and Am 2 / E = 2.4 x 10~ 4 eV 2 /GeV and at which 20' = 0.9tt and 
20" = 1.17T. This point lies in the region of solution A, but very close to the border 
line with solution D: the fact that the condition cos(2#^ — 4:9' m ) = is also fulfilled 
is crucial for having Piv^ — » v s ) = 1. As we have already noticed, for each given 
h ^ 30° there are several total neutrino conversion points at large values of sin 2 29, 
at which the phases 20' and 20" are not necessarily equal to n or to odd multiples of 
7r (see Table 3). Thus, the explanation of the enhancement offered by the indicated 
authors is at best qualitative. 



9 The discussion in the indicated article, being rather qualitative is done for the — > v s transitions 
as the authors presume that cos 20^0 and that the matter term dominates over the energy- 
dependent term in the corresponding system of evolution equations. 
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v e — > v s . The atmospheric v e can undergo small mixing angle v e — > v s transi- 
tions with Aw? from the interval (85). For the Earth-core-crossing neutrinos the 
total neutrino conversion (solution A) takes place at small mixing at Am 2 /E = 
3 x 10" 7 eV 2 /MeV (Table 3), or for E = (3 - 26) GeV if Aw 2 is given by eq. (85). 
Such transitions would lead to a reduction of the rate of the multi-GeV e-like events 
at —1 < cos9 z ^ — 0.8 in the Super-Kamiokande detector. 



5.2 Transitions of Solar Neutrinos 

In the case of the transitions of the Earth-core-crossing solar neutrinos the relevant 
probability is P e 2- The transitions can be generated by v e — v^m or by v e — v s mixing. 
The interference maxima of P e 2, corresponding to solution A & (eq. (0)), P e 2 — 1) 
can take place only in the intervals Am 2 /E 2* (1.9 - 8.2) x 1(T 7 eV 2 /MeV (Fig. 4) 
and Am 2 jE = (0.9 - 3.5) x 10~ 7 eV 2 /MeV, respectively. 

In the case of v e — v^u) mixing, the absolute maximum P e 2 = 1 occurs both at 
small and large mixing angles (Table 4). At small mixing angles it takes place at 
sin 2 20 = 0.038; 0.044; 0.060; 0.077 for h = 0°; 13°; 23°; 30°, respectively. These 
values of sin 2 29 lie outside the region of the MSW SMA solution. The change of P e 2 
in the variables sin 2 29 and Aw? jE along the "ridge" of local maxima leading to the 
"peak" P e2 = 1 at fixed h is illustrated in Table 6 (see also Figs. 1 - 2 in [0] and Fig. 
4b in 0). At sin 2 29 = 7 x 10~ 3 , for instance, we have at the "ridge" P e2 = 0.36; 0.28 
for h = 13°; 23°. These local maxima dominate in P e 2. the local maxima associated 
with the MSW effect in the core (mantle) are by the factor ~ (2.5 — 4.0) (~ (3 — 7)) 
smaller 0. For the values of Am 2 = (4 - 9) x 10" 6 eV 2 of the MSW SMA solution 
region, the dominating local maxima take place at values of the solar neutrino energy 
from the interval E = (5.6 — 13.8) MeV Jl|, to which the Super-Kamiokande, SNO 
and ICARUS experiments are sensitive. 

At large mixing angles, sin 2 29 ^ 0.5, we have P e 2 = 1 at values of sin 2 29 which fall 
in the region of the large mixing angle (LMA) solution of the solar neutrino problem 
(Table 4). For the LMA solution values of Am 2 = (1 - 10) x 10~ 5 eV 2 the absolute 
maxima P e 2 = 1 occur at E ^ 12.5 MeV (see also 0) and typically at even larger 
values of the neutrino energy. Nevertheless, their presence has an influence on the 
transitions even at smaller E and can enhance them somewhat (Fig. 7, see also 0). 

The solar v e can undergo SMA MSW transitions into v s for Am 2 = (3.0 — 8.0) x 
10~ 6 eV 2 (see, e.g., |23|, |10|]). Taking into account the interval of values of Aw? jE 



for which we can have total neutrino conversion, one finds that solution A can be 
realized only for E ^ 9 MeV. 

The implications of the new enhancement effect for the interpretation of the data 
of the solar neutrino experiments have been discussed in [[J and in much greater 
detail in || || || |[|. Let us just mention here that due to this enhancement it would 
be possible to probe at least part of the Aw 2 — sin 2 29 region of the SMA MSW 
—* v ii{t) transition solution of the solar neutrino problem by performing selective 
D-N effect measurements [R |8|, |[|. The Super-Kamiokande collaboration is already 
exploiting successfully these results fTTJ . 



Recently the Super-Kamiokande Collaboration has published for the first time 
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data on the event rate produced only by solar neutrinos which cross the Earth core 
on the way to the detector ||26| : in contrast to their previous night bin N5 data, the 
new Core data are not "contaminated" by contributions from neutrinos which cross 
the Earth mantle but do not cross the Earth core. Due to the new enhancement 
effect, the Super-Kamiokande Core data and the corresponding value of the Core 
D-N asymmetry are particularly useful in performing effective tests the MSW SMA 
v e v h(t) solution, as was suggested in f| and the results of a simplified analysis 
show ||. More specifically, at 2 s.d. a large subregion of the MSW SMA v e — > io r ) 
solution region is incompatible with the D-N effect data (including the Core one), 
while at 3 s.d. the data on the Core D-N asymmetry alone excludes a non-negligible 
subregion of the indicated solution region (see Fig. lc in ||). The results obtained in 
H indicate also that the data on the Core and Night D-N asymmetries can be used 
to perform rather effective tests of the MSW LMA v e — > z/ m(t ) solution as well. 

The consequences of the new enhancement effect for the interpretation of the solar 
neutrino data are less significant in the case of the MSW v e — > v s solution of the solar 
neutrino problem || . Nevertheless, using the Super-Kamiokande data on the Core D- 
N asymmetry |2(J one can probe and constrain also the MSW v e —>■ v s "conservative" 
solution region ||. 

Similarly, the future data on the Core and Night D-N asymmetries in, e.g., the 
charged current event rate in the SNO detector f27j can be used to perform equally 
effective tests of the MSW solutions of the solar neutrino problem M. 



6 Conclusions 

We have investigated the extrema of probabilities of the two-neutrino transitions 
v ii iye) - > v e (^V;t), ^2 — ► ^e? v e — > u s , etc. in a medium of nonperiodic density distri- 
bution, consisting of i) two layers of different constant densities, and ii) three layers 
of constant density with the first and the third layers having identical densities and 
widths which differ from those of the second layer. The first case would correspond, 
e.g., to neutrinos produced in the central region of the Earth and traversing both 
the Earth core and mantle on the way to the Earth surface. The second corresponds 
to, e.g., the Earth-core-crossing solar and atmospheric neutrinos which pass through 
the Earth mantle, the core and the mantle again on the way to the detectors. For 
both media considered we have found that in addition to the local maxima corre- 
sponding to the MSW effect and the NOLR (neutrino oscillation length resonance), 
there exist absolute maxima caused by a new effect of enhancement and correspond- 
ing to a total neutrino conversion, P a p = 1. The conditions for existence and the 
complete set of the new absolute maxima were derived. These conditions differ from 
the conditions of enhancement of the probability of transitions of neutrinos propa- 



gating in a medium with periodically varying density, discussed in [ITH] . The absolute 
maxima associated with the new conditions of total neutrino conversion are absolute 
maxima corresponding to P a p = 1 in any independent variable characterizing the 
neutrino transitions: the neutrino energy, the width of one of the layers, etc. It was 
shown that the new effect of total neutrino conversion takes place, in particular, in 
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the transitions (u e ) — > u e (^ ;T ), ^2 — > v e , v e — > i/ a and P M — > z/ s in the Earth of 
the Earth-core-crossing solar and atmospheric neutrinos. The strong resonance-like 
enhancement of these transitions discussed in Q] is due to the new effect. This 
enhancement was associated in |Q with the NOLR and we show that in certain cases, 
like the v<i — > z/ e = — > z/ e and z/ e — > z/ M ( r ) transitions at small mixing angles, the 
NOLR describes approximately the enhancement. However, in a number of cases the 
NOLR interpretation of the enhancement fails completely. A "catalog" of the most 
relevant absolute maxima corresponding to a total neutrino conversion, was given for 
the transitions indicated above. We have shown that the NOLR and the newly found 
enhancement effect are caused by a maximal constructive interference between the 
amplitudes of the neutrino transitions in the different density layers. Thus, the max- 
ima they produce in the neutrino transition probabilities are of interference nature. 

Our analyses were done under the simplest assumption of existence of two- neutrino 
mixing (with nonzero-mass neutrinos) in vacuum: u e — or v e — v s , or else — v s . 
However, in many cases of practical interest the probabilities of the relevant three- 
(or four-) neutrino mixing transitions in the Earth reduce effectively to two-neutrino 
transition probabilities for which our results are valid. Thus, our results have a wider 
application than just in the case of two- neutrino mixing. 

We have discussed also briefly the phenomenological implications of the new effect 
of total neutrino conversion in the transitions of the solar and atmospheric neutrinos 
traversing the Earth core. The effect leads to a strong resonance-like enhancement 
of the 1/2 — > v e transitions of solar neutrinos if the solar neutrino problem is due to 
MSW small mixing angle v e — > transitions in the Sun. The same effect should be 
operative also in the — > u e [y e — > v^( T )) transitions of the atmospheric neutrinos 
crossing the Earth core if the atmospheric and P M indeed take part in large mixing 
vacuum <-> i> T , «-> v T oscillations with Am 2 ~ (1 — 8) x 10 -3 ) eV 2 , as is 
strongly suggested by the Super-Kamiokande atmospheric neutrino data, and if all 
three flavour neutrinos are mixed in vacuum. The existence of three-flavour-neutrino 
mixing in vacuum is a very natural possibility in view of the present experimental 
evidences for oscillations of the flavour neutrinos. In both cases the new effect of total 
neutrino conversion produces a strong enhancement of the corresponding transition 
probabilities, making the transitions observable even at rather small mixing angles 
(see also |], 0). Actually, the effect may have already manifested itself producing 
at least part of the strong Zenith angle dependence in the Super-Kamiokande multi- 
GeV /i— like data |T|, @, [|. The new effect should also be present in the <-> v s 
transitions of the atmospheric multi-GeV z^'s both at small, intermediate and large 
mixing angles, if the atmospheric neutrinos undergo such transitions. 

As the results obtained in refs. [0, 0, ||, |] and in the present study indicate, 
it is not excluded that some of the current or future high statistics solar and/or 
atmospheric neutrino experiments will be able to observe directly the resonance-like 
enhancement of the neutrino transitions due to the effect of total neutrino conversion. 
Solar neutrino experiments located at lower geographical latitudes than the existing 
ones or those under construction (see the article by J.M. Gelb et al. quoted in [T3|]) 
would be better suited for this purpose. An atmospheric neutrino detector capable 
of reconstructing with relatively good precision, e.g.,, the energy and the direction of 
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the momentum of the incident neutrino in each neutrino-induced event, would allow 
to study in detail the new effect of total neutrino conversion in the oscillations of the 
atmospheric neutrinos crossing the Earth. 
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Table 1: The values of Am 2 /E < 20 x 10~ 7 eV 2 /MeV and of sin 2 29, at which 
the absolute maxima of P efM = P^ e , P es and P^, corresponding to a total neutrino 
conversion, P af3 = 1 (solutions A, eq. ([2"§D), take place for neutrinos traveling from 
the center of the Earth to the Earth surface. The values of 2<j)', 20" (in units of n), 
sin 2 29' m , sin 2 29'^ and sin 2 (26^ — 26' m ) at the absolute maxima are also given. 



P(v^ -> v e ) = P(u e -> u^ t) ) 


^[lO- 7 ^] sin 2 29 


2<j//lT 2<f)"/7T 


sin 2 2C sin 2 2C sin 2 (2^ - 2d' m ) 


6.611 .138 
5.825 .651 
11.736 .649 
18.970 .828 
13.589 .938 


.867 .888 
1.094 1.866 
2.356 2.679 
4.153 4.853 
3.063 3.933 


.434 .603 .999 
.999 .501 .521 
.872 .985 .221 
.935 1.000 .072 
1.000 .885 .117 


P{u e -> u s ) 


3.355 .359 
6.073 .883 
10.181 .904 
13.663 .984 
17.809 .962 


.531 .608 
1.333 1.652 
2.281 2.723 
3.155 3.844 
4.089 4.912 


.774 .862 .590 
.993 .943 .103 
.975 .999 .037 
1.000 .983 .019 
.989 1.000 .012 




4.006 .275 
5.966 .850 
10.441 .849 
13.589 .981 
18.002 .936 


.625 .684 
1.287 1.679 
2.303 2.709 
3.133 3.864 
4.106 4.896 


.612 .745 .870 
.989 .848 .234 
.945 .996 .086 
1.000 .959 .045 
.974 1.000 .027 
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Table 2: The values of Am 2 /E < 20 x 1(T 7 eV 2 /MeV and of sin 2 26, at which the 
absolute maxima of P eA1 = P^ e , corresponding to a total neutrino conversion, P e ^ = 1 
(solution A, eq. (|5"4])), take place for neutrinos crossing the Earth along trajectories 
passing through the Earth core: h = 0°; 13°; 23°; 30°. The values of 20', 20" (in 
units of 7r), sin 2 29' m , sin 2 29'^ and sin 2 (2$^ — 4:9' m ) at the absolute maxima are also 
given. 



h° 


Am' z r-j n-7 eV' z ] 


sin 2 26 


207vr 


20"/7T 


sin 2 29' m 


sin 2 20£ 


sin 2 (2^ - 49'J 


0. 


7.201 


.034 


.923 


.949 


.111 


.614 


1.000 


0. 


4.821 


.154 


.509 


2.334 


.749 


.208 


.296 


0. 


11.256 


.547 


2.174 


4.693 


.795 


.996 


.285 


0. 


3.407 


.568 


.65 


3.612 


.846 


.160 


.881 


0. 


18.113 


.803 


3.933 


9.129 


.923 


1.000 


.710 


0. 


11.292 


.841 


2.425 


6.167 


.987 


.891 


.716 


0. 


19.307 


.903 


4.314 


10.382 


.979 


.987 


.845 


o 


4 587 


948 


1 167 


4 705 


792 


285 


996 


0. 


19.862 


.994 


4.628 


11.738 


.992 


.900 


.981 


0. 


11.373 


.994 


2.702 


7.577 


.954 


.708 


.981 


13. 


7.024 


.039 


.931 


.950 


.131 


.553 


1.000 


13. 


4.537 


.169 


.496 


2.258 


.845 


.179 


.133 


13. 


10.071 


.379 


1.887 


3.180 


.645 


.999 


.068 


13. 


10.954 


.626 


2.272 


4.503 


.869 


.973 


.383 


13. 


3.267 


.654 


.732 


3.456 


.777 


.153 


.966 


13. 


6.736 


.802 


1.477 


4.059 


.997 


.580 


.688 


13. 


15.802 


.874 


3.653 


7.680 


.977 


.972 


.791 


13. 


18.249 


.843 


4.205 


8.603 


.949 


.997 


.762 


13. 


10.741 


.927 


2.529 


5.936 


.998 


.797 


.859 


13. 


18.770 


.979 


4.542 


9.914 


.999 


.922 


.954 


23. 


6.460 


.051 


.910 


.922 


.195 


.387 
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Table 3: The same as in Table 2 for the probabilities P es and Pp$ and for Am 2 /E < 
10 x 10~ 7 eV 2 /MeV. 
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Table 4: The same as in Table 2 for the probability P e2 (solution A & , eq. ([74]) ) in the 
cases of v e — v^u) and v e — v s mixing, and for Am 2 /E < 10 x 10~ 7 eV 2 /MeV. 
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Table 5: Values of P efl = P Me along the "ridge" of local maxima in the variables 
Am 2 /E (in units of 10~ 7 eV 2 /MeV) and sin 2 29 < 0.05 for fixed h, when sin 2 29 de- 
creases from the value at which the absolute maximum of P e ^ = 1 (i.e., total neutrino 
conversion) takes place. The results correspond to Earth-core-crossing neutrinos: 
h = 0°; 13°; 23°. The values of sin 2 26' m , sin 2 20^, sin 2 (20^ - 49'J and 20', 20" (in 
units of 7r), at the local maxima are also given. 
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Table 6: The same as in Table 5 for the probability P e 2 {y e — v^t) mixing). 
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FIGURE CAPTIONS 

Figure 1. The regions of the four different solutions A, B, C and D (eqs. 
([57]) and (T45])) for the maxima of the transition probability P efJL = P^ e in a two-layer 



medium. The two different layers correspond to the core and the mantle of the Earth 



Figure 2. The same as in figure 1 for the case vacuum - Earth mantle. 

Figure 3. The regions of the three different solutions A, C and D (eqs. fl54|), ( |63D and 
(|B~4])) for the maxima of the transition probability P efl = P Me in a three-layer medium. 
The three different layers correspond to the mantle-core-mantle of the Earth. 

Figure 4. The same as in figure 3 for the transition probability P e 2 in the case of 
v e - z^( T ) mixing. 

Figure 5. The same as in figure 3 for the transition probability P^ = P{pn ~^ v a ). 

Figure 6. The probability P efl = P Me for the Earth-center-crossing (atmospheric) 
neutrinos (h = 0°), as a function of sin 2 29 (horizontal axis) and Am 2 /E [10~ 7 eV 2 /MeV] 
(vertical axis). The ten different colors correspond to values of P efl in the intervals: 
0.0 - 0.1 (violet); 0.1 - 0.2 (dark blue); ...; 0.9 - 1.0 (dark red). The points of total 
neutrino conversion (in the dark red regions), P efl = 1, correspond to solution A, eq. 
([54]), for the Earth-core-crossing neutrinos. 

Figure 7. The same as in figure 6 for the probability P e 2 in the case of v e — v^m 
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mixing, and for (solar) neutrinos crossing the Earth (core) along the trajectory with 
h = 13°. The points of total neutrino conversion (in the dark red regions), P e 2 = 1, 
correspond to solution A, eq. (fT4|). 

Figure 8. The same as in figure 6 for the probability P^g and for neutrinos crossing 
the Earth (core) along the trajectory with h = 23°. The points of total neutrino 
conversion (in the dark red regions), P^ g = 1, correspond to solution A, eq. fl54|). 

Figure 9. The same as in figure 6 for the probability P es and for neutrinos crossing 
the Earth (core) along the trajectory with h = 0°. The points of total neutrino 
conversion (in the dark red regions), P es = 1, correspond to solution A, eq. fl54f ). 
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